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Abstract. The Wightman function for a massive free scalar field is studied within 
the light front formulation, while a special attention is paid to its mass dependence. 
The long lasting inconsistency is successfully solved by means of the novel Fourier 
representation for scalar fields. The new interpretation of the light front singularities 
as the high momentum phenomena is presented and adequate regularizations are 
implemented. 



1. Introduction 

The light front (LF) formulation of the relativistic quantum field theory, which starts 
from the pionieering paper by Dirac on forms of the relativistic dynamicspQ, has been 
extensively studied during last twenty years ( for some reviews see j2] , jHj , |1] ) • In 
this approach one introduces the LF coordinates usually as = ± x^)/V2 and 
then picks as the LF temporal evolution parameter, while x~ is treated as the LF 
surface coordinate. The standard procedure, where x^ = t is the evolution parameter, 
refered to as the equal-time (ET) approach, differs substantially from the LF procedure. 
The triviality of LF physical vaccum usually is given as the main advantage of the LF 
approach. 

However when one considers a free field theory then literally no difference should 
appear for all physical quantities defined in either procedure. Especially, the Wightman 
functions should be the same but, as it has been observed already in 1977 by Nakanishi 
and Yabuki |5., it is not the case. They have argued that the ET Wightman function 

^^^^"^^ - (27r)3 J 2uj{pf 



-co 



where Lj{p) = ^rn? + p"^ and x° = t, is a smooth mass- dependent function at the LF 
surface x"*" = 

lim /^et{x) = ^Ki{mJx\). (2) 



They have opposed the LF Wightman function 

oo oo r I 

^ ^0 -oo 

which as long as a;"*" 7^ coincides with - this can be shown be performing explicitly 
all momentum integrations. Then taking the limit x~^ ^ gives the expected result 
(121). The inconcistency arises when one takes x'^ —>■ before doing the momentum 
integrations, since then one obtains 

00 00 

A+^(x+ = 0,x-,xx) = 7^Jy^ J dVe-'[^-^"-^-^'-l. (4) 

^ ^ -00 

The above expression is explicitly mass independent and the integral over p_ variable 
is ill defined. This situation is unexceptable since the value of the Wightman function 
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at the LF surface (x"*" = 0) is crucial for the LF quantization procedure. Further one 
may beheve that the mass dependence and the LF singularity are quite close related 
phenomena. 

The singularities at p_ = are permanent difficulties within the LF formulation and 
different regularizations are used for keeping them under control. However neither the v 
theory of [H] or the DLCQ approach of ^ can fix the above mentioned mass-dependence 
inconsistency. For the free scalar fields the mass appears only via p+ = — then one 
may force that x'^ coordinates of fields should not coincide for the two point Wightman 
function. Since this would stay in a confiict with the LF quantization at the fixed x'^ 
surface, one may allow for the imaginary part of x"*" ^ x"^ — in (jHl) , which effectively 
leads to the formula [3] 



dp_ 







Thus, as long as e > 0, one has both the regularization for small values of p_ and the 
mass dependence. Another solution is even more drastic - one should forget about the 
LF quantization [Sj. Quite recently there was an attempt to solve this mass dependence 
inconsistency within the distribution theory jU], but the conclusion that in the sense of 
distributions one can define only objects which are mass-independent misses the point. 
The partial success of all these approaches can be connected with the LF dogma that 
one has to solve the singularity at p_ — > first, then all other problems will be properly 
solved. We think that one may reverse the logic and start with a different mass- 
dependent Fourier representation of scalar massive field at the LF surface, at least 
for the free field case. 

Our paper is organized as follows. In section |2l we present the standard LF formulation 
the Wightman function for the free scalar massive field. In section IHl we present the 
modified LF formulation. In section |3] we discuss the mass dependence of modified 
LF Wightman function. In section we present different UV regularizations. At 
last, the conclusions and possible further research are presented. The notation and 
technical details of calculations are given in [Appendix A| In [Appendix B we present LF 



canonical formalism for higher derivative Lagrangian. In [Appendix C[ we show how the 
LF singularity may arise within the ET formulation. 

2. Standard light front formulation 

We start our consise presentation of the standard formulation of the light front field 
theory with the free real masive field (^(x), which when quantized on the LF surface 
(x"*" = 0) is represented as the Fourier integral 

(27r)3 io 2k 

where h.c. stands for the Hermitian conjugate and 

= ^J{kl + m^)/2. (7) 



(6) 
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The anihilation and creation operators a(k_|_,fc_ 
commutator 



i^(k_|_,fc_) have a nonvanishing 



a(k_L,fc_),a^(p_L,p 



The Wightman function A+(x) is defined as the vaccum expectation value 

A+(x) = (O|0(x)0(O)|O), 
which due to (jHl) and (jH]) has the Fourier representation 



(27r) 



2fc. 



(8) 



(9) 



(10) 



The integral over the longitudinal momentum must be treated with a high caution 
since for x"*" = (x~ = 0) it diverges logarithmically at the upper (lower) limit of 
integration. Therefore we see that the LF Wightman function is ill defined for the 



space-like separation of points 



< 0, contrary to the ET result 



Thus in order to perform consistently all integrations in (jlO|) we need to keep ^ 0. 
Since for the free field case the Wightman function is known for arbitrary separation 
of points, one may accept this as a kind of technical details which are specific for the 
LF formulation. However from the more general perspective the problem is much more 
serious, since from (jH)) one usually infers the Fourier representation for an interacting 
theory at the LF quantization surface x~^ = 

d^fc I d/c- 



[X 



(2vr) 



2k. 



-ik 



■Xj^^—ik-X 
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Thus we see that any solution of the LF Wightman problem for the free scalar field may 
have a serious consequences for other LF models with interacting scalar fields. Finally 
we would like to present the 4-dimensional integral representation of the LF Wightman 
function 

d2fc. 



[X 



(2vr) 



dA;+dA;_e-''^-^e(fc_)(5(2fc+fc_ - k'i - m"). 



(12) 



which evidently is equivalent to (fTIH) and it explicitly shows that there is no symmetry 
between k± momenta. 



3. Modified light front formulation 

We propose to take another Fourier representation for free real scalar field 

a(k_L, /c_) + h.c. 



[X 



+ 



dk± 


roo 


dA;„ 


(27r)3 




2A;_ 


dA;^ 


roo 


dA;+ 


(27r)3 




2A;+ 



-ifcx-xxg-ifc+x+g ^k+- b{k±,k+) + h.C. 



where the nonvanishing commutators are 

a(k^,fc_),at(p^,p_)] = {2nf2k,6ik,-p^)6''ik^~p^), 
6(k^, A;+),6t(p^,p^)l = {2nY2k_5ik+-p^)S\k^-p^). 



(13) 
(14) 

(15) 
(16) 
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This leads to the modified LF Wightman function 

which can be equivalently rewritten as 

A+(x) = J J dk+dk_e-'''-''eik_ - fXk)S{2k+k_ -kl- m^) + 

+ 1 1 dA;+dfc_e-^'=-"e(fc+ - /ifc)5(2A;+A:_ -k\- m^), (18) 

thus our modification restores the symmetry k^ ^ k^. 

Both integrals in (jl7|) diverge logarithmically in their upper limits when either 
= or x"*" = 0, respectively. This leads us to the conclusion that there are possible 
two UV divergences, contrary to the standard formulation where there is one UV and 
one IR divergency. This new interpretation seems to be more physical since when fields 
are massive no IR divergency should appear. 

When x"*" 7^ 0, we may change the integration variables 

k+ = ^, (19) 

in the second integral and return to the standard result ()1U|) . 

The free field representation (jl3j) can be taken as a basis for the new canonical 
representation (at 



+ 



0) 






dks_ 


POO 


dk_ r 


(27r)3 




2k_ ^ 


dk± 


POO 


dk^ 


(27r)3 




2k+ 



g-ifcx-^xe-^fc-^' a{k^,k_) + h.c. 



e ''=^-^^e b{k^,k+) + h.c. 



+ 



(20) 



for an interacting scalar field. Thus we may compare two LF integral representations 
(|lip and (|20p for their interpretations of modes, which do not depend on x~ . In the 
former one, which is currently commonly accepted, these modes are low momentum 
A;_ = and should be treated as other IR contributions (e.g. by putting the system 
into a finite box in x~ coordinate). In the latter one, which is novel, these modes are 
high momentum (fc+ oo), thus one should treat them in an adequate way (e.g. by 
the Pauli-Villars regularization or by the higher derivative terms). 

We may check the consistency of our modified formula ()17|) by inspecting the LF 
commutator function which is defined as 



A(x) = A+(x) - A+(-x) 
d^/c I /""^ dk_ 



l^-ik-x _ ^-ik-x \ ^-ik±-x± _|_ 



(27r)3 J^(k^) 2A;_ 



Mass dependence of Wightman function 



6 



Now we notice that the k± integrals are convergent, therefore we may easily perform 
these integrations 



d^fc, 

-sgn[x 



- zsgn(a;-) J ^52(/x(fcx)|x-|)e-''=-^- = -^sgn(x-)5^(xx), (22) 

where we have used relations ()1.16|) - ()l.lcj) and we conclude that the commutator 
function has the proper form. 

4. Mass dependence of modified LF Wightman function 

Since the mass dependence (or rather mass independence) of the LF Wightman function 
has lead to serious problems concerning the consistency of LF field theory [S] , thus here 
we will concentrate on the LF surface = for the modified LF Wightman function 

d'^kj_ r'^ dk. 
(2^ A. 2k. 

The simplest criterion, whether our modified LF Wightman effectively depends on m 
or not, is to calculate its derivative with respect to m 
d 



, -A+(x) = / e-^'^'^^ e"^'^^-^^ 

dm 





r d^k± 1 








~ 2"y 


(2vr)3/if 




,x~m 


f d^fc_L f°° 


dA;+ - 


2 J 


' (27r)3 4 





"2 _ 

where we use the simple relation 

-T-f^k = —■ (25) 
dm 

Though our integral ()23p is logaritmically divergent, then its parametric derivative is 
already convergent and the integration over k^ is elementary 

m r d'^k± 1 ^ / i„ _ \ 

+ 2y(27r)3^f 

= --/S^-r^e— (26) 



. {2Txfm? + kl 

Thus quite unexpectedly we find that the final integral does not depend on the LF 
coordinate x~ but only on the transverse coordinates x±. More we see that the mass 
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dependence of the modified LF Wightman function (j23p is by no means triviaL The 
2-dimensional intergal over can be performed in the cyhnder coordinates 

k2 = gsincj), = gcoscj), (27) 

X2 = psinO, X3 = pcos9, (28) 

where the angle integration gives 

''^"d0e'^''™'^(*-'') = 27rJo(^p), (29) 

with Jo{x) being the Bessel function. Then the last radial integral gives 

;^A+ X = - — — / f?df?^^ = -—Ko{mp), 30 
dm \^'^) "'0 + Q'^ 47r^ 

where -^'0(2;) is the modified Bessel function and we have found a well defined mass 

differential equation for the Wightman function. As a quick consistency check for our 

calculations we may take the Wightman function Q and find that it satisfies the above 

mass differential equations, due to the property of the modified Bessel functions 

^[xK^{x)] = -xKo{x). (31) 
ax 

However one may raise the objection, that if one changes the integration variables 

K = (32) 
in the second integral of (j221), then one gets the usual LF representation 

which is evidently mass independent. However, since here this second integral is 
divergent for k^ —* 00, thus such a change of the integration variables is not legitimate. 
Instead, one should first regularize this evident UV divergence, by means of some UV 
regularization and only then, one may change the integration variables. 



5. UV regularizations 

In this section we will define different regularized Wightman function and our starting 
point is the cuttoff expression 

+ 7 (27r)3 4 2A;/ " ^^^^ 

We have subtracted the divergent contribution in a way, which allows to perform easily 
integrations over k± 

A^^(x)= -i/^e-'=---(7 + log(/.,|x-|)+i|sgn(x-)+log^), (35) 



Mass dependence of Wightman function 



8 



where we have used the relations (|l.la|) - (ll.lcp . Further calculations will depend on the 
regularization method that we will choose to get rid of the cutoff parameter A. 

5.1. Pauli-Villars regularization 

The Pauli-Villars regularization is the method of doing perturbative calculations in 
both gauge and Lorentz-invariant way, where one uses the regularized propagators. For 
instance, the massive scalar propagator is 

^'f\^) = / 7^ , , . + E ... . . e---, (36) 
with the following conditions upon the coefficients Ca'- 

N N 

l + ECa = 0, m2 + ^aM2 = 0, etc. (37) 

a=l a=l 

However this regularization is by no means restricted to the pertubation theory, but 
on contrary it can be applied to almost all problems where high momenta divergencies 
appear. For the free scalar field, this method has been applied in the original paper 
[TU] . where Pauli and Villars start with the regularization of singular functions A(x) 
and A'^^)(a;), which are closely related to the Wightman function A"'"(x) 

A+(x) = ^ (A(x) + i^^^\x)) . (38) 

So even in the ET formulation the Pauli-Villars regularization can be applied for the 
free massive scalar field. 

Since the LF Wightman function depends on the cutoff A logarithmically, then the 
Pauli-Villars regularization of this function should have only one auxialiary field (with 
Mi,Ci). Accordingly we define the Pauli-Villars regularization of the LF Wightman 
function as 

where Aik = \l {k\ + Mf )/2. In the current case, only the first condition should be 
taken from (|H7jl 

1 + Ci = 0, (40) 
and further we may use our former results for the integrations over k± given in 



Mass dependence of Wightman function 9 
Since we may introduce the integral representation 

into the integrand, then integrals over the transverse momenta become Gaussians 

A-^(x) = \r'^[ ^e-'=---e-t'=i fe-t-^ - e-f^^H , (43) 
2 JO a J [2'Ky ^ ' 

which lead to 

A-(x) = -i^ ^e-i (e-t-^ - e-t^n = (44) 
^ K,{m\[xi) - ^,^K,{M,Jxi) (45) 



47r2 ri" V 47r2 



As the last step we may push the auxialiary mass to infinity Mi oo and this leads to 

1 TD I 

lim Al'^(a;) = —^KiimJxl), (46) 

\/ _L 

which is the expected result (j2I)- 
5.2. Higher derivative regularization 

Another popular regularization is based on the modified Lagrangians, where one adds 
the higher derivative (HD) terms which effectively change the high momenta behaviour 
of the perturbative propagators. However here, we will use this method for a free field 
theory. The LF canonical procedure for this model is presented in [Appendix B and 



here we will just quote the crucial points. The Wightman function within the HD 
regularization is defined as 

Dl^{x) = (0 |0(x)0(O)| 0) = ^ ((0 |$+(x)$+(0)| 0) + (0 |$_(x)$_(0)| 0)) , (47) 

where we have used ()2.16aj) - ()2.16fe|) . ()2.24aj) . Next we may take ()2.25|) and write 



VI + 4am2 J (27r)3 1^^,+ 2A;+ 2k+ 

(48) 

which is very similar to the Pauli-Villars regularization| ()39p. thus omitting here the 
intermediatery steps, we write down the solution 

D^^{x) = ^ ^ ^ (m+Ki{M+J^) - M_Ki{M^J^S\ . (49) 



I One should not be surprised that the HD resuhs are very similar to those from the Pauli-Villars 
regularization, since one may prove that, for some models, these methods are equaivalent [Tl|. 
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where we write after (I2.22al) - (l2.226j) and (jT^ 

= i-(l - Vl + 4am2) ^ m^, (50a) 



M! = — (1 + Vl + 4am2) oo. (506) 
2a 

Therefore we may take the hmit a — 0, removing the HD regularization, 

777 I 

hm Df^ix) = =KUmJx^, ), (51) 



which again coincides with 

6. Conclusions and further prospects 

In this paper we have proposed the modification of the LF Wightman function, which 
properly solves the problem of mass dependence. Our modification starts with the 
Fourier representation for the free massive scalar field operator and perfectly agrees 
with the canonical commutation relations. Further we find that there is no longer any 
IR singularity problem in this LF field theory, but rather a new UV divergency arises. 
We point out that the distribution sgn(x~) in the commutation function (j^^ is gener- 
ated by the high momentum behaviour of the LF Wightman function. We may further 
speculate that for an interacting theory one should take another Fourier representation 
(r20|l . which may lead to new description of the LF systems. 

The analogous analysis for the higher spin fields, specially fermions and gauge fields, 
will be given elsewhere, since these cases are technically more complicated and here we 
have decided to present our conjuctere for the simplest case of quantum field theory. 
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Appendix A. Notation and useful integrals 

We use the natural units c = h = 1. Our LF notation starts with the definitions of null 
components for the coordinates = (x° ± x^)/-\/2, while the transverse components 
). The similar definitions are taken for any 4- vectors. The LF surface 
coordinates are denoted as a; = (a;~,a:;*). The partial derivatives are taken with respect 
to contravariant coordinates, thus we have d+ = d/dx^,d^ = d/dx~,di = d/dx\ The 

metric tensor has non vanishing components g^^ = g ^ = l,gij = The scalar 

product of 4-vectors is a ■ 6 = a_|_6_ + a_6_ — aj6j, while for the LF surface components 
we have a -b = a^h_ — aibt. 
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In the main text, the following integrals are helpful 

/ — cos/e+ / — (cos /c — 1) = — 7 — log/i 
Ja k Jo k 



and 



— sinkx = -sgn(a;) [tt — 2Si(/i|x| 
^ k 1 

/ — — sin — X = sgn(x)Si(/i|x|)., 
Ju k k 



:i.la) 



[1.1b) 



:i.lc) 



where Si{x) is the integral sine function. In [Appendix C we use the following integral 



dggMgp) 



2m 



l-<5 



which is valid for 6 > 1/4. Also we use another type of integral 



da 



m 



^i7r5/2sgn(it) 



a 



1-5 



\u\ 



;i.2) 



;i.3) 



which is valid for < i?e 5 < 1. 

Appendix B. Simplified LF canonical formalism for higher derivative 
Lagrangian 



Let us describe the free scalar field with the higher order derivative and start with the 
Lagrangian density 



2 " 2 " 2 
which leads to the Euler-Lagrange equation of motion 

, 2" 



0, 



(2.1) 
(2.2) 



where we use d"^ = d^d'^ and a > is the regularization parameter. 

For the sake of canonical quantization procedure, we prefer to extend the field 
content of our system by adding another scalar field x which leads to the equivalent 
Lagrangian density 



a 2 



_02 -^-d^, 

2 ^ 2a ^ 
with the system of Euler-Lagrange equations of motion 



a 



0. 



In the LF formulation, we explicitly we write 

1 Tn'^ x^ 

C'^ = d+<pd-<p - -{di<pf - d+<pd-x - d-cpd+x - di(pdiX - —<P^ - — • 



(2.3) 

(2.4) 
(2.5) 

(2.6) 
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The component of the energy-momentum tensor is defined as 



5L' 



+ ^^d+x = {d- 



— 2c?_0(9_x- 



(2.7) 



(59+0 5d+x 

where L'^ = J d^xC'^. Here we will use the recently proposed simplified LF canonical 
quantization procedure [13^, which is based on the trivial equations 



id_(j){x) = (f){x), P'^{x^) 
id^xix) = [x(x),P+(x+) 
with the translation generator 



(2.8) 
(2.9) 



P+(x^ 



d-^xT^ 



d-^x 



{d^(j){x) - d^x{x))^ - (d^xix))^ ■ 

From the diagonal structure of P'^, we immediatelly recoqnize fields 
independent LF canonical fields, with nonvanishing LF commutators 

2 [0(x) - x{x), d.(f){y) - d-x{y)]x+=y+ = ^^H^ - v), 
'^[x{x),d^x{y)l^+=y+ = -iS%x-y), 

or equivalently 

2[<p{x),d^x{y)l+=y+ = -^S\x-y), 
'^[x{x),d^x{y)L+=y+ = -iS^{x-y), 
[<P{x),d^<P{y)l^^y^ =0. 

Now we will look for the independent modes for our system and suppose that they 
are given as the linear combinations 

$+ = c+0 + x, (2.16a) 

<|._ = c_0 + x, (2.16&) 

while their equations of motions are 

(2a+a_ - + Ml) $± = 0. (2.17) 

These expressions are compatible with the equations of motion ()2.4|1 and ()2.5j) 

X 



(2.10) 

X and X as the 

(2.11) 
(2.12) 

(2.13) 
(2.14) 
(2.15) 



(29+a_ - Ai)0 
(29+a_ - A^)x 



a 
X 

a 



+ m 



provided the following relations are satisfied 

,2 



Ml 



m 
1 + c± 
a 



with the consistency condition 



0. 



If we choose the regularization parameter to satisfy the following inequalities 
1 



4^2 



< a < 0, 



(2.18a) 
(2.186) 

(2.19a) 
(2.196) 

(2.20) 
(2.21) 
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then 



-1 + VI +4am2 

c+ = < 0, (2.22a) 



-1 - Vl + ^OiW? 

< 0, (2.226) 



2 

and we end up with two nontachyonic modes, since we have 

2 

Tfl 

Ml = -— > 0. (2.23) 

Now we may find the LF commutators for the independent modes 

$+(a;+,x),<l>_(a;+,|/)] = - i {c+ + c_ + 1) 6^{x - y) = 0, (2.24a) 
<l>+(a;+,x),$+(x+,?/)l = - i (2c+ + 1) 6^{x - y) = -Wl + Aam?5'^{x - y), (2.246) 



<l>_(x^,a;), <l>_(x^,x) 



i (2c_ + 1) 5^(x - y) = Vl + 4:am^%x - y), (2.24c) 



which indicate that $_ and $+ are the positive and negative metric fields, respectively. 
Thus the nonvanishing LF Wightman functions for the independent modes can be 
defined, accordingly to (fTTjl . as 

A+(x,a) = (0|$„(x)$,(0)|0) = 

7 (27r)'^ J^^fc, 2/c_ 

+ aVl + 4am2 / ^e-'^+^'e-^-^-e-^-", (2.25) 

J [271) Jfi^^a 2/c_|_ 



with (a = ±) and Hk± = V(M| + kl)/2. 

Appendix C. ET momentum integrals in cylinder coordinates 

The LF singularities that we have discussed in the main part of this paper, can be found 
also within the ET quantum field theory. In order to encounter them one has to perform 
the momentum integrations in (0) in the cylinder coordinates and in d space dimensions 
[Ej. Thus we take coordinates (pi, g = ^Jpl + . . . + p^, i^d-i) and the volume element is 

d'^p = g'^-'dg dp, dn,_,, (3.1) 

where Qd-i is the solid angle in d-1 transverse directions. Then the angular integral 
over transverse directions is 

/2 \ ('^-3)/2 

dfi,_ie'^^-^^ = 271 i—j Jid~3)/2{QP), (3.2) 



where p = y x^ + ■ ■ ■ + x% 

Now we may concentrate on the integration over variable pi, 



I{t,X^) = J ^e-^i-^t-P^^') = - J d^g-iV^^^(tcosh^-xismh„)^ ^3 3^ 
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where we have introduced rj variable through the parameterization 

Pi = \J w? + smhrj. 
The result depends on the value of [T^ 

oo 

— CO 

oo 

- j (\.Tje~^^ m?+Q'^\fs^ cosh Tj _ 
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(3.4) 



for < 0, 



-oo 
oo 



lY^^Hl"\^m^ + Q^^) for s2 > 0, 



for t = , 



(3 



-oo 
oo 



for t = —x^ 



where a = 1 for t > 0, and cr = 2 for t < and evidently at the LF surfaces = this 
integral is ill-defined. 

Let us consider the singularity at x° = x^ and we introduce the modified integral 



l\t,x' 



drje 



iy^ m^+(5^(t cosh?;— X"*- sinh rj) —rj —Srj 



(3.6) 



with < 5 < 1. Thus we find 



I\t,t) = - / dr^e 



r(5); 





1 ^isgn(f)<57r/2 



2|t|5(m2 + ^2)5/2 

The remaining integration over q is simple 



(3.7) 



oo 

/ df?f?^+V,(f?p)^^ 



+ ^2)5/2 



l-<5/2 



r(5/2) 



Ks/2-i-e{mp), (3.8) 



for 5 > e + 1/2, where e = (c/ - 3)/2. 

Now we would like to remove our regularization parameter 6 by taking the limit 
5^0, but this is allowed only if e < —1/2 and this is the reason why we have started 
with the dimensional regularization d < 3. Thus in d dimensions we have the ET 
Wightman function 



. (27r)i+^ /mV+'"^/' r((5) e 



5^0 {2ttY \ p 
m f m ^ ^ 



T{6/2) 2\t 



isgn{t)5n/2 



Ki+,{mp). 



(3.9) 



47r2p V27rp^ 

At last, we may go to physical three space dimensions and obtain for any space-like 
distances ^2 < 



A^(x; m^) 



m 



Att'^P 



Ki{mp). 



(3.10) 
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Though our final result has the commonly known form (j2)), we stress that we have 
encountered in our ET calculations similar problems to those, which appear within the 
LF formulation. The cylinder coordinates clearly single out one space direction, thus 
two light fronts = appear naturally. In contrast, when one uses the spherical 
coordinates all space direction are equivalent and no LF singularity arises. 
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